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Transient Uncoupling Induces Synchronization 
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Finding conditions that support synchronization is a fertile and active area of research with 
applications across multiple disciplines. Here we present and analyze a scheme for synchronizing 
chaotic dynamical systems by transiently uncoupling them. Specifically, systems coupled only in 
a fraction of their state space may synchronize even if fully coupled they do not. Although, for 
many standard systems, coupling strengths need to be bounded to ensure synchrony, transient 
uncoupling removes this bound and thus enables synchronization in an infinite range of effective 
coupling strengths. The presented coupling scheme thus opens up the possibility to induce synchrony 
in (biological or technical) systems whose parameters are fixed and cannot be modified continuously. 


Synchronization is one of the most prevalent collec¬ 
tive phenomena in coupled dynamical systems [1]. Syn¬ 
chronization and related consensus phenomena have been 
frequently found in biological, ecological, physical, engi¬ 
neering and social systems such as in predator-prey dy¬ 
namics, the spread of epidemics, the migration of large 
populations, systems of self-driven particles and systems 
of social or technical dynamics [2-12]. For chaotic sys¬ 
tems, synchronization typically emerges only within a 
specific range of coupling strengths and is impossible oth¬ 
erwise [1, 13-15]. 

In this letter, we propose and analyze a way of 
inducing synchronization between coupled chaotic 
oscillators by transient uncoupling: If the system is in 
a certain predefined subset of its state space, coupling 
is active; otherwise it is inactive. We systematically 
study the dependence of successful synchronization on 
the fraction of state space where coupling is active. 
Synchronization may emerge even for systems that 
coupled continuously in time (i.e., standard coupling) 
do not synchronize. Furthermore, the system may 
synchronize for an infinite range of coupling strengths, 
although this is often not possible for ordinarily coupled 
chaotic systems. A systematic numerical analysis 
reveals how transverse stability properties vary across 
the attractor with the location of active coupling, not 
only between more or less stable synchrony but all the 
way from stability to instability for the same system. 
This demonstrates that transient uncoupling modifies 
the collective dynamics in a non-trivial way. These 
results may find applications in inducing synchrony 
in systems whose local coupling parameters cannot be 
continuously varied with ease but only switched on or off. 
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Standard coupling. To start, consider a system of 
two unidirectionally coupled chaotic oscillators 


dxi 

dt 

dX2 

dt 


F(xi), 

F(x 2) -b aC X (xi - X 2 ) , 


( 1 ) 

( 2 ) 


where Xi (t) , X 2 (t) G denote the states of the driv¬ 
ing and driven unit, respectively, C is a square cou¬ 
pling matrix, and a the coupling constant which deter¬ 
mines the overall strength of coupling ]13]. As an ex¬ 
plicit example throughout this Letter we consider iden¬ 
tical x-coupled Rossler oscillators defined by F(x) = 
(—(j/-b z), cc-b ay, 6-b z(x — c))^ ]16] and C € 
where Cy = 1 for i = j = 1 and = 0 otherwise. Fur¬ 
ther, a = b = 0.2, c = 5.7 and we take x^ =: Zi)~^ 

as a convenient notation. Other chaotic systems exhibit 
qualitatively the same phenomena as those presented be¬ 
low ]17]. 

Depending on the coupling strength a, such systems 
do or do not synchronize towards xi(t) = X 2 (t) =: xs(t). 
In particular, like many other coupled chaotic systems, 
Rossler oscillators are known to typically synchronize for 
intermediate coupling strengths a, but not if coupled too 
strongly or too weakly (Fig. la-c). 

These qualitative synchronization properties depend 
on the (‘transverse’) dynamics of the difference xj^ = 
xi — X 2 . A Taylor expansion to first order in the (xj_)j 
yields 


XX = F(xi) - F(x 2 ) - aC X (xi - X2) 

« ]J(xs(t)) - aCjxx (3) 

where J(x) = 9xF(x) is the local Jacobian of F. For the 
system (3) to relax to x±{t) —>■ 0, its maximum transverse 
Lyapunov exponent 
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needs to be negative ]1]. Fig. Id illustrates a 

function of the coupling constant a. This clearly links in 
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FIG. 1. Synchronization depends on coupling strength. Tra¬ 
jectories of the driving (solid line) and driven (dashed line) 
unit of two coupled chaotic oscillators for (a) a = 0.05, (b) 
a = 1.5, and (c) a = 5 as indicated in panel (d). (d) The 
maximum transverse Lyapunov exponent A,(iax indicates syn¬ 
chronization for intermediate coupling only. 


a quantitative way the coupling strength and the qual¬ 
itative changes in collective dynamics observed before 
(Fig. la-c). 


Transient uncoupling. We now introduce transient 
uncoupling via a factor 


Xa(x2) 


1 for X 2 G A ; 
0 for X 2 ^ A , 


(5) 


in the coupling term, 

^ = F(x 2)-baxA(x2)C X (xi - X 2 ). (6) 

Here H C is a subset of the driven unit’s state space 
where coupling is active. The two units are thus effec¬ 
tively coupled only within a subset x H of their com¬ 
mon state space. For A = the units are ordinarily 
coupled continuously in time. 

Practically relevant subsets A are defined by clipping 
a region of state space along the direction of a particular 
coordinate axis, 

Aa = {x2 G : |(x2)i - (x2)i| < A}, (7) 

where X 2 is a suitable point and the subscript 'V refers 
to the first coordinate of X 2 and x^. Thus, coupling is 
only active within a column of width 2A centered around 
(x 2 )i. Here, (x 2 )i = 1.2 was chosen as the center of 
the attractor in a;-direction. An example realization for 
Rossler oscillators is illustrated in Fig. 2. 

Such transient uncoupling modifies the collective dy¬ 
namics of the coupled system in a non-trivial way (Fig. 3 
and 4). Specifically, for a fixed coupling strength a, for 
which standard coupling would not lead to synchroniza¬ 
tion, clipping in an intermediate interval A induces syn¬ 
chronization. Obviously, for A —>■ 0 the units become 


I -2A- 1 

X* - A X* X* + A 



FIG. 2. Transient uncoupling through state space clipping. 
The dynamics of two synchronized chaotic oscillators in the 
x-y plane with x* = (x 2 )i = 1.20, A = 4.16, and a = 7.0 
(driving: solid curve; driven: dashed curve). Coupling is only 
active in the interval x-i G \x* — A,®* -|- A] (shaded in gray). 


completely uncoupled and cannot synchronize. Similarly, 
for no clipping A —0/2 (where O is the width of the 
attractor along the clipping direction) we reobtain the 
original system with standard coupling which does not 
synchronize. For intermediate clipping, however, we find 
stable synchronization. As the clipping fraction A be¬ 
comes just one additional parameter of the system we 
expect the Lyapunov exponent to vary continuously with 
respect to A. An analysis of the transverse Lyapunov ex¬ 
ponent as a function of the clipping fraction A' = 2A/n 
confirms this (Fig. 3). 

Intriguingly, we find that for a fixed clipping interval 
A, the dependence on the coupling strength a is changed 
not only quantitatively but also qualitatively (compare 
Fig. 4 to Fig. 1). In particular, for intermediate tran¬ 
sient uncoupling (intermediate values of A) synchrony 
emerges in an infinite range of coupling strengths a, 
thus in particular for arbitrarily large coupling (Fig. 4d). 
This is in contrast to many chaotic oscillators which, 
when ordinarily coupled, exhibit an upper bound above 
which synchronization fails [13]. In fact we explicitly 
checked that the same phenomenon also emerges in 
Rossler oscillators with other parameters and in pairs of 
coupled Lorenz and coupled Chen oscillators as well as 
for larger networks [17]. 

Optimal uncoupling. We now analyze direction de¬ 
pendencies of transient uncoupling. Interestingly, the 
range of coupling strengths a for which the system syn¬ 
chronizes increases when the clipping fraction decreases 
from A' = 1, as Fig. 5a illustrates. Moreover, the range of 
clipping fractions for which synchronization emerges de¬ 
pends on the exact direction in state space along which 
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FIG. 3. Synchronization induced by transient uncoupling. 
Maximum transverse Lyapunov exponent for two transiently 
uncoupled chaotic oscillators (parameters see text) for a = 5 
and clipping with (x 2 )i = 1.20. Synchronization emerges for 
moderate clipping, i.e., intermediate values of A', although 
not without clipping (A' = 1). 




FIG. 5. Extended synchronization range by transient un¬ 
coupling and optimal clipping, (a-c) Depending on the cou¬ 
pling strength a and the percentage (A' = 2A/fI) of the 
state space where coupling is active, the system may or may 
not synchronize. The dark area marks the parameters where 
the synchronized state is stable, i.e., AjJiax < 0. Clipping 
is (a) in x direction {x* = (x 2 )i = 1.2), (b) in y direction 
{y* = (x 2)2 = —1.5), and (c) in the direction y « 3.la: 
{9* = O.dyr). Accordingly, the direction of clipping can be 
optimized to achieve the largest possible clipping range, (d) 
Effectiveness S{9) (Eq. 8) of clipping along the direction 6 for 
fixed a = 10. 


the fraction of clipping 

s{e)= [\{f,e)df (8) 

Jo 

for which the system synchronizes when a is fixed. Here, 
we have measured the angle 9 counterclockwise from the 
a:-axis and have defined the synchrony indicator 


FIG. 4. Transient uncoupling induces synchronization in an 
infinite range of coupling strengths. Trajectories of the driv¬ 
ing (solid line) and driven (dashed line) units for (a) a = 0.05, 
(b) a = 1.5, and (c) a = 5, the same as in Fig. 1. The clipping 
is given by Eq. 5 and 7 with (x 2 )i = 1.20 and A = 4.16 as 
in Fig. 2. (d) Maximum transverse Lyapunov exponent 
as a function of the coupling strength a, note the logarithmic 
scale. The grey line shows A^ax for normal, undipped cou¬ 
pling. With transient uncoupling synchronization is stable for 
arbitrarily large coupling strengths. 


clipping is applied. For instance, clipping along the x- 
axis seems more synchronizing in this sense than clipping 
along the y-axis (compare Fig. 5a to Fig. 5b). Oblique di¬ 
rections exhibit even broader ranges of clipping fractions 
where synchrony emerges (Fig. 5c). 

In fact, certain directions of clipping are optimal. Due 
to the shape of the attractor, excursions of trajectories 
that substantially vary z are rare compared to those that 
vary the other two coordinates. Thus clipping is desirable 
in the x-y plane. To quantify the effectiveness of clipping 
depending on its direction in the x-y plane, we measure 


s{f.e) 


1 for ^max < 0 ; 

0 for ^max > 0 , 


and the temporal clipping fraction 

/= lim 1; f XA(x2(t))dt, 


(9) 


( 10 ) 


such that larger values of S{9) indicate that synchroniza¬ 
tion emerges in a larger range of clipping fractions. 

The curve S{9) has two local maxima (Fig. 5d), 
indicating two locally optimal clipping directions, one of 
which is globally optimal (at 9* « 0.47r). Why is there 
such a complicated dependence on direction? 


Transverse stability depends on uncoupling lo¬ 
cation. We can better understand the synchronization 
of the system by characterizing the stability for a family 
of clipping functions x(x) = Xt(x) where 

A = := {x 2 e : |x 2 - x^l < r (x^)} . (11) 

Coupling is thus active if and only if X 2 is in a sphere 
of radius r (x^) around x^. We sample the center points 
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FIG. 6. Multiple switches between stability and instability 
depending on the coupling location. The color of the points 
indicates the maximum transverse Lyapunov exponent of the 
system with clipping to sets A = j4x*,r (Eq. 11), indicating 
stable synchronization (light) and no synchronization (dark) 
depending on the coupling location on the attractor. Param¬ 
eters are a = 5 and / = 0.05. 

randomly from the attractor (i.e., the invariant measure) 
of the uncoupled system and choose the size r (x^) such 
that the coupling is active during a fraction / of the time. 

The results show that the impact of the uncoupling 
strongly depends on the position where clipping is ap¬ 
plied. In particular, at identical system parameters, syn¬ 
chrony can be either stable or unstable, depending on 
where the coupling is active (Fig. 6). This holds even 
though the coupling is active for the same fraction / 
of time. The attractor regions of positive and negative 
transverse Lyapunov exponents alternate depending on 
the direction from the origin. As these different regions 
of stability and instability each occur two times on the 27r 
phase cycle (circulating the origin) and at roughly equal 
phase distance, this explains the two maxima (and the 
two minima) of the curve S{9) found above (Fig. 5d). 
This heterogeneous dependence on the exact location in¬ 
dicates that transient uncoupling, despite being repre¬ 
sented by a linear reduction of the coupling term, mod¬ 
ifies the collective dynamics of the system in a strongly 
nonlinear way. As a consequence, the clipping sets A 
need to be determined individually for each given system 
to be synchronized. 


any upper bound on the coupling strengths that enable 
synchronization. As a natural extension, it would be 
challenging to explore how systems capable of weaker 
forms of collective dynamical coordination, such as phase 
synchronization, lag synchronization or generalized syn¬ 
chronization would behave if transiently uncoupled [1]. 
Additionally, our scheme may extend synchronization 
regimes not only in continuous-time systems (described 
by differential equations, discussed throughout the Let¬ 
ter), but also for chaotic maps and systems temporally 
switching between different continuous dynamics, cf. 
e.g., [22, 23]. 

Stability properties of chaotic systems are known to 
vary locally with the system’s state as quantified by 
the local Lyapunov exponent [24-26]. For transverse 
systems, studied above, local stability depends on the 
direction of the difference vector Xj^ = Xi — X 2 . For 
small coupling strengths the direction of this vector in 
the uncoupled transverse system accurately indicates 
the regions of state space where coupling will be most 
effective. However, when the coupling is stronger or 
active in an extended region of state space the trajec¬ 
tories are more strongly modified by the coupling. In 
particular, whether coupling at one point is effective or 
not in general depends non-linearly on the coupling in 
the rest of state space. Optimizing the regions of active 
coupling in this respect might enhance synchronizability 
even further. 

As experimental chaotic systems often exhibit intrin¬ 
sically fixed or at least restricted internal and coupling 
settings, the question emerges how to synchronize them. 
Transient uncoupling by state space clipping may help 
to induce synchronization for a wider range of coupling 
strengths, with potential applications to chaotic lasers, 
electric and electronic circuits, communication systems 
and chaos based cryptography [27-37]. 
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